Plasma diagnostics using digital holographic interferometry 
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The advances in Charge Coupled Devices in one hand and the high resolution measurements of 
holographic technique on the other hand, we have adopted the method of digital real-time holo- 
graphic interferometry for the diagnostics of high density plasma. The measured values of plasma 
electron density agree with the measurements from other techniques. 



In recent years, digital holographic interferometry and 
shadowgraphy are emerging as a versatile tool for non 
destructive testing and diagnostics purposes in materials 
engineering Similar techniques are also adopted for 
diagnosing plasma and to measure various plasma param- 
eters using pulsed lasers. However a plasma diagnostics 
method that can reliably and rapidly measure fluctua- 
tions with both spatial and time resolution can lead to in- 
creased understanding and eventual control of turbulent 
transport. Additionally, reliable diagnostics that mea- 
sure the plasma profiles are critical to moving towards 
an actual prototype magnetic fusion reactor. Conven- 
tional methods for plasma diagnostics include [1| Lang- 
muir probe Q, microwave, laser interferometry and 
Thompson scattering [1, Out of these, laser interfer- 
ometry has the advantage as it is versatile and gives ac- 
curate results. Since the advent of optical lasers, inter- 
ferometry has been widely used to study different types 
of plasmas j6l| . We have made an attempt to successfully 
adopt the digital holographic interferometry for plasma 
diagnostics. To authors knowledge, this is the first time 
digital holographic interferometry is employed for plasma 
discharge tubes in real-time. This open up the way for 
the diagnosis of arc plasma, DC discharge plasma, etc. 
The advantages of this technique over other conventional 
spectroscopic methods are its high resolution measure- 
ments of the order of few tens of nanometer, real-time 
measurements and cost effectiveness. Further this meth- 
ods could give realistic parameters of plasma like electron 
density and electron temperature 

The present technique provides information in the form 
of two-dimensional maps of the electron density, without 
the need of extensive modeling. In a laser interferome- 
try various types of interferometer can be used such as 
Michelson interferometer, Mach-Zhender interferometer 
(MZI) etc. We employed a Mach-Zhender interferom- 
eter. Chord integrated phase information is utilized to 
obtain the radial profile of electron density by using Abel 
inversion technique with appropriate software code writ- 
ten using MathCAD^M and LabVIEW"^". 

We adopt the fact that the geometrical and optical 



lengths traveled by light through plasma are different, 
since the refractive index of plasma is proportional to 
the density of free electrons. When a laser beam passes 
through a cylindrical tube containing circular symmetric 
plasma density (rig), the light suffers a change in optical 
path length. Light passing through different chords of the 
tube as shown in Fig. 1, are phase shifted by different 
amount. The change in phase is estimated from 0,3 




FIG. 1: Schematic of light propagation inside plasma tube 
of radius R. The solid horizontal arrows parallel to x axis 
represents laser beams. The geometrical path length traveled 



by light inside plasma tube at ith chord is 2xi = 2^J~i^ 
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where yi is the horizontal distance of ith chord from ori- 
gin, Aq, e, Arieiyi), c, and eo are the wavelength 
of laser light, charge of an electron, change in electron 
plasma density at ith chord, speed of light, mass of elec- 
tron and absolute permittivity, respectively. Free elec- 
trons produce negative phase shift because the effective 
index of refraction is less than 1. Hence, the change in 
refractive index will be less than 1, but positive. 

The experimental setup is shown in Fig. 1. A 30 mW 
polarized He-Ne laser (Melles Griot) having central wave- 
length at 632.8 nm passes through a spatial filter assem- 
bly and a beam expander (BE). A cylindrical lens is used 
to convert the circular symmetric beam into a line beam. 
The line helps us to measure the optical properties along 
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a single cross section of the plasma tube only. Also, the 
plasma tube, if mounted on a linear stage may be use- 
ful to scan the plasma tube completely. The expanded 
beam is divided using a 50:50 non-polarizing cubic beam 
splitter (BSl). A part of the beam is allowed to inci- 
dent on a plane mirror (hereafter, we call it reference 
beam) while the other part is reflected by another mir- 
ror (object beam). In this experimental setup, we have 
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FIG. 2: Schematic of experimental setup used for digital holo- 
graphic interferometry. NDF - Neutral density filter, BE - 
beam expander, CBSl, 2 - cubic beam splitters, Ml, 2 - mir- 
rors, PTl, 2 - high pressure Hg-plasma tubes, CCD - charge 
coupled device. 

used high pressure mercury lamps (Phillips) as a plasma 
source. Two identical lamps are used in both the arms 
of the Mach-Zhender setup. One lamp acts as a compen- 
sator while the other lamp is connected to power source. 
Both the beams interfere inside a non-polarizing cubic 
beam splitter (BS2). Interference fringes are monitored 
using a charge coupled device (CCD) (Apogee, Model- 
LISAA-M). A large area convex lens assembly (Pentax) 
is used to collect the beam into the CCD. The whole 
system is mounted on a vibration isolation table (Melles 
Griot). Interference fringes are viewed from one arm of 
BS2 using a screen. In the absence of any plasma, the in- 
terference fringes (or chord integrated intensity) recorded 
in the CCD can be expressed as 



^{V^) =h+l2 + 2^hhcos{(j>). 



(2) 



Here, Ii and I2 are the intensity of the laser in the two 
arms of the interferometer while <j){t) is the phase differ- 
ence between the two beams. Since, the two beams are 
passing through a path which is identical but intention- 
ally at a small angle so as to obtain straight line fringes. 
In the presence of plasma eq. may be rewritten as 



/(?/„ t) = /i + /2 + 2^/hhcos[^{t) + A(j){y„t)]. (3) 

Simulated results obtained using eq. ([3]) and exper- 
imental observation of chord integrated intensity are 
shown as contour plots in Figs. 3a and 3b, respectively. 
Due to the limitation of the experimental setup we could 



obtain only a partial image from the CCD. A typical in- 
terference signal obtained using CCD is shown as Fig. 
3c. The fringes continuously obtained with time are dis- 
played in Fig. 3b. While obtaining the simulated results, 
the values of constants are obtained from the present 
experimental conditions. The experimental data is ob- 
tained from CCD, while the simulated results are gener- 
ated from a software code written using Lab VIEW. For 
times t = P to 50 sec no current is applied to PTl. Since 
no temporal change in phase is occurring, the fringe pat- 
tern remains almost same as evident from Figs. 3a and 
3b. At t — 50sec, constant current is applied to the 
plasma source. In the presence of plasma the interference 
pattern are modulated with time. The number of fringes 
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FIG. 3: (a) Contour plot of simulated results chord integrated 
intensity with time and spatial positions, (b) Contour plot of 
experimental results of chord integrated intensity with time, 
(c) Chord integrated intensity at t = llSsec is shown with 
pixel positions. In all the figures, the position of chord is 
represented as yi. 

increases due to large change in plasma density, but sat- 
urates after the build-up (rising) time of plasma. The 
unknown parameter Ane{yi) for different chords is ob- 
tained after performing many iterations of the simulated 
results for minimum standard deviation. Best profile of 
Ane{yi) which matches gives minimum standard devia- 
tion with the experimental profile is used for calculation. 

The measurement of phase change indirectly gives the 
magnitude of electron density. Figure 4, shows the tem- 
poral change in plasma density. When the plasma tube is 
switched ON, the electron density raises rapidly and sat- 
urates. However, it cools slowly after switching it OFF. 



3 





t— — I — ■ — I — ■ — I — ■ — h- 



04.25 
03.30 



02.55 



1.275 



20 40 

time (sec) 



80 



The chord integrated phase changes with respect to 
axial distance (j/i) are converted radial distance by us- 
ing Abel inversion method and are shown in Figure 5. 
The axis on the left shows the numerical values of phase 
change while the axis on the right side is the estimated 
values of change in electron density after Abel inversion. 
These values agree with the recent observations using 
spectroscopic method fTo(. 



FIG. 4: Contour plot of change in phase with axial position 
and time. The plasma tube is switched ON at 18th second 
and switched off at 35th second and allowed cool till 90sec. 
Maximum phase change of Air is recorded. The bottom curve 
is data at the center (r = 0mm). 




-6-3 3 

radial / axial distance (mm) 

FIG. 5: Experimental measurement of phase change (ob- 
tained from chord integrated intensity) with axial distance 
and electron density (estimated from emissivity profile ob- 
tained after Abel inversion process) with radial distance. 



To conclude, we used the method of holographic inter- 
ferometry to estimate the plasma electron density. The 
measured value agrees with the standard values of high 
pressure mercury plasma density. The present work has 
the potential of measuring, imaging the plasma density 
variations in real-time. This method can also be used 
for obtaining the 3D profile of electron density in a DC 
plasma. 
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